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Abstract. - We study the transition at T = from a ferromagnetic insulating to a ferromag- 
netic metallic phase in manganites as a function of hole doping using an effective low-energy 
model Hamiltonian proposed by us recently. The model incorporates the quantum nature of 
the dynamic Jahn-Teller(JT) phonons strongly coupled to orbitally degenerate electrons as well 
as strong Coulomb correlation effects and leads naturally to the coexistence of localized (JT 
polaronic) and band-like electronic states. We study the insulator-metal transition as a func- 
tion of doping as well as of the correlation strength U and JT gain in energy Ejt, and find, for 
realistic values of parameters, a ground state phase diagram in agreement with experiments. 
We also discuss how several other features of manganites as well as differences in behaviour 
among manganites can be understood in terms of our model. 



The colossal magneto-resistance(CMR) exhibited by manganites {Rei-xAxMnOs , Re = 
La, Nd, Pr etc. and A = Sr,Ca, Ba etc.) for a range of hole doping x around x ~ 0.3 and 
near the Curie temperature Tc, where they undergo a transition from a low temperature fer- 
romagnetic metallic phase to a high temperature paramagnetic insulating phase, has led to 
a great deal of interest ^12] in these systems, which also show a variety of other interesting 
phenomena such as charge and orbital ordering and incipient phase separation. The interplay 
of orbital degeneracy of the itinerant Cg electrons of Mn, their coupling to lattice degrees of 
freedom, especially to degeneracy removing Jahn- Teller (JT) phonons, strong Coulomb cor- 
relation effects and related Hund's rule coupling (between the eg electrons and the t2g core 
spins of Mn) are believed to be responsible for these phenomena, but achieving a detailed 
theoretical understanding has been a major challenge. 

We have recently proposed a new effective low-energy Hamiltonian starting from 

the two qualitatively different coexisting vibron states at the each site of the lattice, one 
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consisting of localized JT polarons ;6^, which we label and the other, which we label 6, 
dispersive and forming a broad band as outlined below. In the presence of strong correlation 
U and ferromagnetic Hund's rule coupling Jh on site between Cg and t2g spins this leads 
O0j to a consistent description of many known features of these systems such as the finite 
temperature insulator metal transition (IMT) near Tc for a range of x, CMR, the existence 
of a ferromagnetic insulating phase at low doping, good metallic behaviour of electron doped 
systems, etc., and a heuristic understanding of several other features, such as the isotope effect 
on Tc and the two-phase coexistence seen by a variety of experimental probes over a range of 
X and T In this paper, we present in detail our theory for the ground state behaviour 

of manganites in the ferromagnetic state as a function of doping x (for x < x^o where charge 
ordering occurs) in terms of the Falicov-Kimball model (FKM) |S] involving the £ and b states. 

The un-doped compound, eg. LaMnO^ is an Mn — O bond (JT) distorted but structurally 
ordered, A-type anti-ferromagnetic insulator with ferromagnetic order of the t2g spins in plane 
(defined by the JT distortion pattern) and anti-ferromagnetic order perpendicular to it. On 
hole doping, say with Sr, anti-ferromagnetism disappears at x = Xci = .08 and the ground 
state is a ferromagnetic insulator till x = Xc = .16, beyond which it is a ferromagnetic metal. 
For LaCa, Xd = 0.1 and Xc ~ 0.18 whereas for NdSr, Xci = 0.18 and Xc is not accurately 
known. The occurrence of the fully ferromagnetic insulating phases and the IMT at such 
large values of x cannot be understood in a model for manganites with only eg electron 
double exchange caused by Jh 01121 • For, since doping generates holes (unoccupied sites) in 
the Cg band, one expects a metal especially when the t2g spin alignment is ferromagnetic, 
since the Cg electrons can then move without hindrance. This is true even for large U. 

We outline below (see ref. PI for details) how the FKM 8 describing the correlated £ and b 
states arises in the context of manganites at T = from a conventional lattice model with two 
Cg orbitals per site, large electron JT phonon coupling g and perfect spin alignment resulting 
from a new ferromagnetic 'virtual double exchange ' coupling Jp ■ We treat it using dynamical 
mean field theory (DMFT) [10], and show that the b band has a reduced effective width, 2D, 
roughly equal to y/x times its bare width 2Do for large U. Hence, below a critical doping Xc the 
b band bottom lies above the localized JT polaronic £ levels of energy —Ejt, which are then 
the only ones occupied, leading to an insulator. For larger x, the 6-band bottom lies below 
the £ level, so the system is metallic. The ferromagnetism is largely due to Jf (with a small 
contribution from conventional double exchange in the metallic phase). Thus the observed 
T = pattern of phases and phase transitions follows naturally from our picture I I Ij . We 
discuss the ground state properties of the FKM as a function of the hole doping (x) and 
the model parameters, namely U, Dq and Ejt, and show that for realistic values of these 
parameters, the calculated phase boundary is in good agreement with experimental trends. 
We conclude by discussing some other implications of our work in regard to experiments. 
A detailed discussion of the finite temperature properties of our model, including the ferro 
metal- para insulator transition and the CMR, is presented elsewhere [Hj- 

We start with a model with two degenerate Cg orbitals per site and strong degeneracy 
breaking electron JT phonon coupling g. For large g, at each site there will be one vibron 
solution , labelled £, with its energy reduced for single occupancy due to a large JT distortion 
by an amount Ejt — {g"^ /2K) ~ leV ^21 where K is the force constant of the JT phonon 
mode. The orthogonal vibron solution [SI, labelled b , is not JT distorted, and hence has no 
gain in JT energy. Inter-site hybridization of the £ states is reduced by the phonon overlap 
or Huang-Rhys factor ry = exp{—EjT/2huJo), which for manganites is ~ (1/200) since 
{EjT/huJo) ~ 10. As a first approximation, we neglect this altogether and treat the £ states as 
site localized. Inter-site hybridization amongst the b electrons is not suppressed. They hence 
form a broad band. In the presence of large U the b states have their largest amplitudes on 
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hole sites [x) as they are strongly repelled from the polaronic sites t. In the presence of large 
Jh the spins of both £ and h electrons are aligned parallel to the t2g core spins. 

An additional consequence 3 of the existence of JT distorted, localized I states is that 
virtual, adiabatic hopping processes involving them in the presence of large U and Jh (where 
an £ electron at site i quickly hops to an empty neighbour j and back, with an intermediate 
state energy cost of 2Ejt due to the unrelaxed lattice distortion at i) give rise to a new, major, 
doping dependent ferromagnetic nearest neighbour exchange coupling Jp between the t2g core 
spins of order t^{l — x)x/{2EjtS'^). In the ground state for x > Xci when this interaction 
dominates the anti-ferromagnetic super-exchange, the t2g are fully ferromagnetically polarized; 
then so also are the eg (both £ and b) spins due to the large Jh{» t). Hence the spin degrees 
of freedom are frozen. 

Thus we are led to effectively spin-less localized t and mobile b electrons with a strong 
local Coulomb repulsion U . The relevant effective Hamiltonian is just the Falicov-Kimball 
model (FKM) , given by 

Heff = -J2 ^^M^^ - i^j^ + E - A* E ^l^^ + E ""^'^^ (1) 

(ij) i i i 

Here b} and £| create the band and localized polaronic states described above at site i, and 
tij are effective, orbitally averaged ^1] inter-site hopping amplitudes for the b states. The 
chemical potential fj, is determined by the doping- dependent constraint that the total number 
of Cg electrons is (1 — x) per site. The relevant parameter regime of H^ff we are interested in 
for manganites corresponds to i ^ 0.2 eV U very large (~ 5 — SeF ^21) and Ejt in the 
range 0.5- 1.0 HSHig. 

While there are no known techniques for exactly solving H^ff, a dynamical mean field 
theory (DMFT) treatment of it can be carried out exactly jJUj. In this approximation, which is 
exact in infinite dimensions, and quite accurate for three dimensions |1U|. the lattice problem is 
mapped to a single site problem embedded in a self consistent effective medium or electron bath 
that represents all the other sites of the lattice. We assume that the system is homogeneous 
The b electron self energy Ey (a;) due to the interactions U is site local i.e. Sij(w) = 
6ijT,((jj), and is determined from the single site or local effective action 

Seff ^ - / dTdT'b\T)g-\T-T')b{T')+(3{-EjT- ^i)nt + Un, dTUbir) (2) 
Jo Jo Jo 

Here &^(t) and &(t') are fluctuating fermionic Grassmann fields, n£ = 1 or corresponding to 
the £ state being occupied or not, and ^/(t) is the bare on site local propagator in the presence 
of the effective medium. The local partition function, obtained by summing exp(~Seff) 
over all the configurations, is expressible as Ziocai = Zq + Zi. Here Zo,Zi arc constrained 
partition functions corresponding to = 0, 1 respectively, and are calculable using standard 
techniques TU' in terms of g{uj~^), the Fourier coefficients of CJ(t) analytically continued to real 
frequencies. Explicitly, Zq = exp{Pao{T)) ; Zi = exp{Pai{T)) x exp [—/?(—£' jt — /i)]; where 

am{T) = TT~^ J dujnp{uj ~ ^) Im {£n{Q~'^{uj+) ~ U6,ni)} , (3) 

n^(uj) being the Fermi function [1 -I- exp{l3uj)] ^ . The local single particle Green's function 
G(uj^) is given by 

G(w+) = ~{bb+)s^„ = wag{iv+) + wi{g-\iv+) - (4) 
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with wi = Zi/{Zq + Zi) — fii and wq — {1 — wi) being the annealed probabihties for the £ 
state being occupied and empty respectively. The condition that the average Cg occupancy 
ni + fih = — x) determines the chemical potential fi. A re-normalized or effective l electron 
energy can be defined by writing wi = n,p{el — whence e| = —Ejt + cio{T) — ai{T). 
Clearly, at T — 0, fi is necessarily pinned to e| as long as fii is non zero. 

In the DMFT 10 Eq. 01 relating G(w+) and G{(^^) for our model has to be supplemented 
by two other equations: namely the Dyson equation, G~^{uj~^) = Q~^{uj~^) — and 
the self consistency relation which expresses the local Green's function in terms of the lattice 
Green's function, as G([j+) = / deDoie) / {tu'^ + fj. — e — I1{uj~^)) , where £>o(e) is bare density of 
states (DOS) for the b-band. These coupled equations for Q, G and S have to be solved self- 
consistently, and typically numerically, to obtain all the quantities of direct physical interest. 
The self consistency relation becomes algebraic, considerably simplifying such calculations, 
for a semicircular DOS, i.e., for -Do(e) = {'^/'^Dq) y/ — where Dq is the half bandwidth. 
This DOS, exact for the Bethe lattice in infinite dimensions ^UIj is fairly accurate for our 
model in 3d, especially for trends and magnitude estimates |lti) . Hence we confine ourselves 
to the semicircular DOS results in this paper. In this case, G{z) = 2/{z + \/ — D^) where 
z = (w+ + M ~ ^('^^)) is complex. Using this result, Eq. 0] the Dyson equation, and the 
equations for wq^w\ we have numerically solved the DMFT equations self-consistently for a 
wide range of values for the correlation J7, JT polaronic energy Ejt^ bare bandwidth 2_Do and 
doping X. The results are discussed below. A broad perspective of the trends in these can be 
obtained in the simpler limit when [/ — > oo. Hence we discuss these first. 

In the U ^ oo limit, the DMFT equations stated above (for the semicircular DOS) can 
be solved analytically. The local Green's functions have the simple form 

G(w+) = 2wo/ {u;+ + fi+ ^{iJ++fiy-D^'^ = wog{LJ+) 

with D = ^/woDq. The local spectral function or re-normalized DOS of the 6-band is simply 
p(w) = 2a/ D'^ — (uj + fj,)^ / (ttDq^), i.e., once again of the semi-circular form, but with a re- 
duced effective bandwidth 2D, and reduced weight wq. At T = 0, fib, wq — I — fi£ and can 
hence be evaluated from (numerically solving) the equations: 

e| = -Ejt + ififib/wo); 

m = (wo ~x)^ (wo/tt) [sin-^i^i/D] + {7t/2) + {^i/ D)^/T^~ii^/D)^} + D). 

In addition, when ng ^ 0, there is the pinning condition e| = /i. These equations have the 
self consistent solution fib = 0, e| = —Ejt — /i, for Ejt > D = y/woDo — y/xDo, since 
Wq — X ior fib = 0. Thus we have the analytic result that for x < (Ejt/Dq)^ the effective 
half bandwidth D = y/xDo, the localized £ levels lie lower than the b band bottom, only the 
former are occupied and the system is an insulator. The T = electrical gap between the 
occupied £ levels and the unoccupied b band bottom in this ferro-insulator phase is given by 
A = Ejt — D = Ejt — ^/xDo- The critical doping for the T = ferro-insulator to ferro- metal 
transition, determined by the vanishing of A, is thus Xc — {Ejt/DqY. As x increases beyond 
Xc, the system becomes a ferromagnetic metal, with e| the re-normalized £ level lying above 
the b band bottom, fib increases with x ,till at some value Xc2, = (1 — Xc2) so that fi^ = 0. 
Beyond Xc2, only the b band is occupied. 

Thus our theory leads naturally to an insulating (ferromagnetic) state for x < Xc, a ferro- 
metallic state with coexisting band b and localized £ electrons for Xc < x < Xc2, and a metal 
with only b electrons and bandwidth 2Do for x > Xc2- Indeed, for x close to 1 ("electron 
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doped limit") there are surprisingly successful calculations J7] of (magnetic) ground states 
based on a model of independent broad band tight binding electrons moving in appropriate 
magnetic superstructures of t2g spins which have AF super-exchange interactions, completely 
ignoring JT interactions. Our theory provides a rationalization for this. 

Detailed zero temperature results from our theory are shown in Figs. l(a)-(c) and in Fig. 
2. We choose Ejt = 0.5 eV and Dq = leV. Fig. 1(a) shows the variation of the re-normalized 
£ level and the b band edge positions {U = oo, full Hne, U = 5eV, a realistic value, dotted hne) 
with doping x. The effective b bandwidth 2D becomes very small (0 for U = oo) as a; ^ 0. 
The physical reason is that, as mentioned before, the b electrons reside mostly on the hole 
sites (fraction x), being strongly repelled from those occupied by £ polarons (with repulsion 
energy U >> 2D). At x = Xc {— 0.25 for U = oo, and very nearly this value for U = 5eV), 
the band bottom crosses the £ level. The system is metallic for x > Xc, with both n;, and 
being nonzero, as shown in Fig. 1(b). We note that in the 'metallic' regime Xc < x < 0.5, nf,, 
the average band occupancy, is small, eg. at x = 0.4, n;, ~ .08. For the set of parameters 
chosen, the £ level empties out completely for x — Xc2 — 0.72, and beyond this only band 
states are occupied. Double exchange and anti-ferromagnetic super-exchange describe the 
magnetic behaviour of the system in this regime. In Fig. 1(c), the effective band or electrical 
gap is shown as a function of x. The gap vanishes smoothly and rises as a; — *■ to 

a value Ejt, and not U ; this should be the electrical gap seen at any finite x no matter how 
small, or in a T 7^ experiment. Fig. 2 shows how the T = insulator metal boundary (xc) 
shifts as a function of Dq and U. Xc increases as Dg decreases, and as U increases. We note 
that all physical quantities for U = 5eV are close to those for ?7 = oo. 

We now briefly compare our results with observed material trends and experimental num- 
bers where available. One of our predictions is that of a ferromagnetic but insulating ground 
state (FI) for x < Xc- The former arises from the new virtual double-exchange coupling Jp 
described earlier; the state is insulating because the effective b half bandwidth D < Ejt 
for X < Xc- In contrast, the ferromagnetic state is necessarily metallic if it arises solely due 
to double exchange Experimentally, all doped manganites have an insulating, fully po- 
larized ferromagnetic state, for Xci < x < Xc ', eg. for LaCa, Xci = 0.10 and Xc = 0.18. 
In our calculations, Xd — (because we have neglected the small, orbital order depen- 
dent, anti- ferromagnetic exchange important at small x jll)'). while the critical doping Xc 
for the insulator-metal transition at T = depends on material parameters roughly as 
Xc = {Ejt/Dq)'^ (cf. Fig. 2). This prediction can not be directly compared with experi- 
ment since the systematics of Ejt and Dq are not precisely known. It is however believed |18) 
that the bare half-bandwidth Dq decreases in the sequence LaSr, LaCa, NdSr (and PrCa), 
because of cation size and its effect primarily on the Mn-O-Mn bond angle and via this on the 
nearest neighbour hopping, while Ejt does not change much. The observed Xc for this se- 
quence has values 0.16, 0.18, and ~ 0.20, the trend being consistent with our prediction. The 
puzzle as to why some manganites (eg. PrCa) have only insulating ground states unlike the 
above can be understood within our theory in terms of the characteristic values of {Ejt/ Dq) 
appropriate to the materials (eg., from Fig. 2, for Ejt ~ .5eV, the ferro-insulator extends 
up to a; = .5 if Dq - 0.7 eV). 

The electrical activation energy in the Ferro-insulator state goes as Ag// = DQ{y/x^— y/x) 
for large U in our theory. This dependence cannot be compared with experiment, since there 
are no measurements of activation energies close to T = 0; the only experimental results we 
are aware of ^H] are for T > Tc, in the paramagnetic phase. The corresponding gaps are not 
expected to go to zero at Xc. However, this high T(> Tc) activation energy does decrease as 
expected from our theory. 




Fig. 1 - (a) Variation of b band edges and the effective £ level e^, in units of Do, with doping x. 
(b) The number n; (of localized JT polaronic £ electrons) and fib ( of band electrons), per site, as 
a function of x. (c) The T = insulating gap A = {Ejt ~ D), as a function of x. In all the cases 
EjT ~ 0.5 eV and Do — 1 eV. The full line corresponds to [/ = cx) and long dashed line to f/ = 5 eV. 
Fig. 2 - The critical concentration Xc for the T = insulator metal transition as a function of Do the 
bare b half-bandwidth, in eV. The full line corresponds to ?7 = oo and line with circles to U — 5 eV . 
The inset shows Xc as a function of U/Do for Ejt ~ 0.5 eV and Do = 1 eV. The point corresponding 
to U = 5 eV is marked by a filled circle. 

An additional consequence of our model is that in the ferromagnetic metallic ground state, 
the concentration of mobile (6) electrons is rather small (Fig. 1(b)), and not {1 — x), the total 
number of Cg electrons per site. This is exactly the inference from the small Drude weight, 
i.e. the integrated optical conductivity, which is a direct measure of the effective number of 
carriers. For example, Okimoto et. al. |^ find Ueff — 0.06 for LaSr at a: = 0.3. Our results 
for fif, quoted above are very much in this range of smallness. A related consequence of our 
theory, arising from the fact that the large majority of the eg electrons are in polaronic £ 
states even in the metal, is the persistence of local polaronic distortions well into the metallic 
phase, for which there is considerable experimental evidence [H]. Thus our results provide a 
natural explanation for several unusual ground state properties of manganites. 

In summary, we have presented here a new coexisting polaron/broad band electron model 
for manganites, which revives the Falicov-Kimball model in a new, hitherto unexpected and 
unexplored, setting. We have completely solved the model in this new context within the 
DMFT, and shown that this leads to a physical explanation and a quantitative theory of 
many characteristic and hitherto puzzling ground state properties of doped manganites in the 
doping regime 0.1 < x < 0.5. 
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